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ABSTRACT 

Massive black hole binaries, with masses in the range 10 3 — 10 8 M©, are expected to be the most 
powerful sources of gravitational radiation at mHz frequencies, and hence are among the primary targets 
for the planned Laser Interferometer Space Antenna (LISA). We extend and refine our previous analysis 
(Sesana et al. 2004), detailing the gravitational wave signal expected from a cosmological population 
of massive black hole binaries. As done in our previous paper, we follow the merger history of dark 
matter halos, the dynamics of the massive black holes they host, and their growth via gas accretion and 
binary coalescences in a ACDM cosmology. Stellar dynamical processes dominates the orbital evolution 
of black hole binaries at large separations, while gravitational wave emission takes over at small radii, 
causing the final coalescence of the pairs. We show that the GW signal from this population, in a 3 year 
LISA observation, will be resolved into ~ 90 discrete events with S/N > 5, among which ~ 35 will be 
observed above threshold until coalescence. These "merging events" involve relatively massive binaries, 
M ~ 10 5 M Q , in the redshift range 2<z<6. The remaining ~ 55 events come from higher redshift, less 
massive binaries (M ~ 5 x 10 3 M Q at z>6) and, although their S/N integrated over the duration of 
the observation can be substantial, the final coalescence phase is at too high frequency to be directly 
observable by space-based interferometers such as LISA. LISA will be able to detect a fraction >90% 
of all the coalescences of massive black hole binaries occurring at z<5. The residual confusion noise 
from unresolved massive black hole binaries is expected to be at least an order of magnitude below the 
estimated stochastic noise. 

Subject headings: black hole physics - cosmology: theory - early universe - general relativity - 
gravitational waves 



1. INTRODUCTION 

Gravitational radiation, described as a tensor pertur- 
bation to the metric travelling at the speed of light, is 
a natural consequence of Einstein's general relativity. It 
has been recognised (e.g., Thorne 1987) that black hole 
binaries are among the most important sources of gravi- 
tational waves (GW), both for ground based interferome- 
ters such as LIGO (Abramovici et al. 1992) and VIRGO 
(Bradaschia et al. 1990), and for the planned Laser Inter- 
ferometer Space Antenna (LISA, Bender et al. 1994). 

Interferometers operate as all-sky monitors, and the 
data streams collect the contributions from a large num- 
ber of sources belonging to different cosmic populations. A 
precise determination of stochastic GW backgrounds from 
different classes of astrophysical objects is therefore cru- 
cial to interpret the data. While a GW background may 
provide information on the number density, redshift evo- 
lution, and mass function of the emitting population, con- 
fusion noises add to the instrumental noise limiting the 
possibility of detecting other class of objects. Moreover, 
to optimize the subtraction of resolved sources from the 
data stream, it is important to have a detailed description 
of the expected rate, duration, amplitude, and waveforms 
of events. 

LISA will operate in the frequency range 0.01 mHz - 1 
Hz, where GW emission from a cosmological population 
of massive black hole binaries (MBHBs) is expected to be 
important (Haehnclt 1994). Today, massive black holes 



(MBHs) are ubiquitous in the nuclei of nearby galaxies 
(see, e.g., Magorrian et al. 1998). If MBHs were also 
common in the past, and if their host galaxies experi- 
ence multiple mergers during their lifetime, as dictated by 
popular cold dark matter hierarchical cosmologies, then 
MBHBs will inevitably form in large numbers during cos- 
mic history. The formation and evolution of MBHs has 
been investigated recently by several groups (e.g., Menou, 
Haiman & Narayanan 2001; Volonteri et al. 2003), and 
the expected GW signal from inspiraling MBH binaries 
has been first discussed by Rajagopal & Romani (1995), 
and recently by Jaffe & Backer (2003), Wyithe & Loeb 
(2003), Sesana et al. (2004, hereafter Paper I), and Enoki 
et al. (2004). 

In Paper I we computed the GW background from MB- 
HBs and the number of coalescences observable by LISA 
in a 3-year mission, adopting the scenario for the assem- 
bly and growth of MBHs proposed by Volonteri et al. 
(2003a, b). In such model, "seed" holes are placed within 
rare high-density regions (minihalos) above the cosmologi- 
cal Jeans and cooling mass at redshift 20. Their evolution 
is followed through Monte Carlo realizations of the halo 
merger hierarchy combined with semi-analytical descrip- 
tions of the main dynamical processes, such as dynamical 
friction against the dark matter background, the shrinking 
of MBH binaries via three-body interactions, their coales- 
cence driven by the emission of gravitational waves, and 
the recoil associated with the non-zero net linear momen- 
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turn carried away by GWs in the coalescence of two un- 
equal mass black holes (the "gravitational rocket" ) . Major 
halo mergers lead to MBH fueling and trigger quasar activ- 
ity. In this paper we use the same model to provide a more 
detailed characterization of the GW signal from inspiral- 
ing MBHBs. Their contribution to the LISA data stream 
is twofold: unresolved sources will give origin to confusion 
noise to be compared to instrumental noise and other as- 
trophysical stochastic backgrounds (e.g. from white dwarf 
binaries, Farmer & Phinney 2003), while resolved inspi- 
raling binaries will probe gravity in extreme conditions 
(e.g. , Vecchio 2004) . Confusion noise and resolved sources 
should provide different cosmological information. The 
former, produced by a large number of unresolved MB- 
HBs, will trace light MBHBs at very high redshift, placing 
constraints on black hole formation scenarios prior to the 
reionization epoch; the latter will be a formidable tool to 
follow the cosmic evolution of MBHs and the formation 
and dynamics of MBH binaries following galaxy mergers. 

The plan is as follows. In § 2 we review the basics of the 
detection of GW from MBHBs, defining observable quan- 
tities such as the characteristic strain amplitude, signal- 
to-noise ratio, and source detection rate. In § 3 we briefly 
summarize our scenario for the cosmological evolution of 
galaxy halos and associated holes. In § 4 we present con- 
fusion noise levels and source number counts. Finally, in 
§5 we discuss our results. 

2. GRAVITATIONAL WAVE SIGNALS 
2.1. Bursts and periodic events 

Following Thorne (1996), an interferometer can be char- 
acterized by two different sensitivity curves, depending on 
the type of signal one expects to detect, i.e. a "burst" or a 
"periodic" GW source. A burst, a short-lived signal whose 
waveform can be utterly complicated, can be described in 
terms of a characteristic strain amplitude h c at the ob- 
served frequency f c ~ 1/ At s , where At s is the duration of 
the signal (Thorne 1987). The spread of the power spec- 
trum around f c will be A/ ~ f c . At the other extreme, a 
perfectly periodic source emits, for the entire duration of 
the observation, at a fixed frequency /. The power spec- 
trum will be peaked at /, with a spread A/ ~ f/N, where 
N is the number of wave cycles clipped into the observa- 
tion. In this respect, a burst can be thought as a single 
complete waveform with f — f c - In the case of a peri- 
odic signal, the interferometer sensitivity is increased by 
the fact that, across the observing interval r, the signal is 
repeated fr times. 

The sensitivity to bursts (/is) and to periodic signals 
(hp) are related by: 



In Figure 1 the two curves hs and hp are com- 
pared for an assumed 3-year LISA observation. The 
curves are obtained combining the LISA single-arm 
Michelson sensitivity curve (taken from the URL 
www.srl.caltech.edu/~shaiie/sensitivity) with the recent 
analysis of the LISA instrumental noise below 10~ 4 Hz 
(Bender 2003, extended from 3 x 10~ 6 Hz to 1 x 10~ 6 Hz 
with a constant slope). 



Consider now a periodic signal of finite duration, with 
strain amplitude h. The total energy carried by the wave 
will be proportional to the number of wave cycles n spent 
at that particular frequency. The quantity to be compared 
with fig is then the "characteristic" strain h c = hy/n. 
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Fig. 1. — LISA single-arm Michelson sensitivity curve to bursts 
(thick solid line) and periodic signals (thick dashed line) in a 3- 
year mission. Data are obtained from 

www.srl.clatech.edu/~shane/sensitivity, and Bender (2003). The 
strain h (dashed lines) and characteristic strain h c (solid lines) for 
a MBHB with M 2 = O.lAfj = 1O 5 M at z = 1 (upper lines), and 
M2 = Mi = 10 s Mq at z = 7 (lower lines) are also shown. 

Note that for a periodic signal at frequency / lasting 
for a time interval longer than the observation time r, we 
have simply n — fr. Then, the signal-to-noise ratio S/N 
increases by the same factor one would obtain comparing 
h to hp in equation (1). The former approach, i.e. com- 
paring h c to hs rather than h to hp, is more general, as 
it allows us to characterize the S/N not only for perfectly 
periodic signals (n = fr), or for bursts (n = 1), but also 
for events in which the emitted frequency shifts to increas- 
ingly larger values during the spiral-in phase of the binary 
system. In the latter case, n = 11(f) represents the num- 
ber of cycles spent in a frequency interval A/ ~ / around 
frequency /, and hence h c is the strain in a logarithmic 
frequency interval (Flanagan & Hughes 1998). Typically, 
the timescale for frequency shift is long compared to the 
wave period, and short compared to the duration of the 
observation. Only close to the innermost stable circular 
orbit (ISCO), the GW frequency changes at a rate com- 
parable to the frequency itself (n ~ 1 and hence h c ~ h). 
In Figure 1 we also show h and h c for two representative 
binary systems. One should note that the true observable 
GW signal is, for f > n/r (the "knee" frequency in the 
h c curves), lower than h, as for these high frequencies the 
source is not monochromatic over the observation time. 
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Such simply consideration naturally leads to define h c and 
h B . 

2.2. Characteristic strain 

Consider now a binary system at comoving distance 
r(z). The strain amplitude (sky-and-polarisation aver- 
aged) at the rest-frame frequency f r is 



8tt^3 G 5 ^M 5 / 3 2/3 

' ~ 10V2 C 4 r (z) - tr ' 



(2) 



where = (M 1 M 2 ) 3/5 / (Mi + M 2 ) 1/5 is the "chirp mass" 
of the binary, and all the other symbols have their stan- 
dard meaning. The strain is averaged over a wave period. 
The rest-frame energy flux (energy per unit area per unit 
time) associated to the GW is 



dE _irc 3 
dAdt ~ 1~G 



(3) 



As discussed above, the important quantity to consider is 
the number of cycles spent in a frequency interval A/ ~ / 
around a given frequency /. Assuming that the backreac- 
tion from GW emission dominates the orbital decay of a 
binary, during the spiral-in phase one can write 



n * fr/jr = 



■fr 



5/3 



(4) 



967r 8 / 3 G 5 / 3 7W 5 / 3 
where we have used the rest-frame frequency shift rate 

Note that n can be computed either in the rest or in the 
observer frame. The characteristic strain in an observation 
of (observed) duration r is then 



h c = h^/n - 



and 



1 



3 l/ 27r 2/3 c 3/2 r(2 ) 



— 1/6 



n > fr, 



n<fr, (6) 



(7) 



where / = f r /(l + z) is the observed frequency. Using 
Parseval theorem, it is easy to see that h c is related to 
the Fourier transform of the strain h, as h\ = 2f 2 h 2 (f r ), 
where h is defined over the positive frequency axis. The 
specific energy per unit area is then 



dE 



--h 2 

4G c ' 



dAdf, 

and, from equation (6), we obtain 



dE Tr 2 / 3 



_ = 1 (-2/3 .,5/3 f-l/Z 

df r 3 ^ ™ h • 



(8) 



(9) 



Note that dE/df r cx f r 1/3 , while (eq. 3) dE/dt cx # 



10/3 



2.3. Signal-to-noise ratio 

In an operating interferometer, any stochastic signal will 
add up (in quadrature) to hs to form the effective rms 
noise of the instrument, h rms . An inspiraling binary is 
then detected if the signal-to-noise ratio integrated over 
the observation is larger than the assumed threshold for 
detection, where the integrated S/N is given by 



S/N Af 



r/+A/ 



din/' 



hc(ft) 



*(/') 



(10) 



Here, / is the (observed) frequency emitted at the starting 
time t = of the observation, and A/ is the (observed) 
frequency shift in a time r starting from /. The latter is 
implicitly given by 

rf+Af df 



If r ' 



(ii) 



where df/f = (1 + z)df r /f r . The frequency at the ISCO 
is, strictly speaking, defined only in the test particle limit 
M2 <§C M\. In the general case, various estimate of the 
transition point from in-spiral to plunge exist, and differ 
by a factor of 3 at most (e.g., Kidder, Will & Wiseman 
1993; Cook 1994). Such uncertainties do not affect our re- 
sults in any manner, so we use, for the observed frequency 
at the ISCO, the conventional Kcplcrian defintion: 



/iSCO 



1 



6 3 / 2 ttG (Mi + M 2 ) 



(1 + *)- 



(12) 



Replacing / + A/ with /isco in equation (11) gives tisco, 
the time needed to span the frequency interval [/, /isco]> 
to be compared a 'priori to r. In the case r > tisco? 
we then set / + A/ = /isco m equation (10). In Figure 2 
we plot h c for different MBHBs at different redshifts, com- 
pared to the LISA h Tms (see § 4.1) multiplied by a factor of 
5, assuming a 3-year observation. If h c > 5/i rms , then the 
signal has, approximately, an integrated S/N > 5. This is 
for illustrative purposes only, as the actual S/N must be 
integrated over the observing period using equation (10). 

At frequencies higher than the "knee" , the time spent 
around a given frequency is less than 3 years, and h c cx 
j-1/6 The signal shifts toward higher frequency during 
the observation, and reaches the ISCO and the coales- 
cence phase in most cases. The lowest curve represents 
a low mass, high redshift equal mass binary. As we shall 
see below, these sources are common in our hierarchical 
model for MBH assembly. In terms of their detectabil- 
ity by LISA, they represent a somewhat different class 
of events. Contrary to the case of more massive bina- 
ries present at lower z, the final coalescence phase of light 
binaries lies at too high frequecies, well below the LISA 
threshold. 

For frequencies much below the knee, the characteris- 
tic strain is proportional to / 7 ^ 6 , as the timescale for fre- 
quency shift is longer than 3 years. The signal amplitude 
is then limited by the observation time, not by the intrinsic 
properties of the source. The source will be observed as a 
"stationary source" , a quasi-monochromatic wave for the 
whole duration of the observation. An increase in the ob- 
servation time will result in a shift of the knee toward lower 
frequencies. The time needed for the sources to reach the 
ISCO starting from the knee frequency is, approximativcly, 
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the observing time. Figure 2 shows that very few station- 
ary sources above threshold should be expected anyway. 




-6 -4 -2 

log observed frequency [Hz] 



Fig. 2. — Characteristic strain h c for MBHBs with different 
masses and redshifts. From top to bottom, the first three curves 
refer to systems with log(Mi/M0) = 7,6,5, respectively, and 
M2 = O.lMi. The solid, long-dashed, and short-dashed lines as- 
sumes the binary at z = 1,3,5, respectively. A 3-year observation 
is considered. The lowest solid curve assumes an equal mass binary 
Mi = M2 = 10 3 Mq at z = 7. The small diamonds on each curve 
mark, from left to right, the observed frequency at 1 year, 1 month 
and 1 day before coalescence. The thick curve is LISA 5/i rms (see 
§ 4.1), approximatively the threshold for detection with S/N > 5. 



2.4. Coalescence rate and number counts 

Given a coalescence rate R, using the frequency shift 
rate / (eq. 5), we can solve for the mean number of indi- 
vidual binaries resolved during an observing period r. We 
begin considering that a MBHB spans, during its lifetime, 
a finite frequency range, / m in < / < /isco, where the 
lower limit is set to the observed frequency at the hard- 
ening radius (Quinlan 1996). Then, from continuity, the 
number of individual observable MBHBs can be computed 
as 

/•/isco Jl 

N T = R % + Rt. (13) 

The first term is simply the integrated density of sources 
in the frequency domain, and does not depend on r. It 
is the number of sources caught in a snapshot of the en- 
tire sky. The second term is the number of new binaries 
born (at frequency / mm ) during the observation time r, 
and must be equal to the number of coalescences within 
the same period. 

The general argument above does not consider that real 
detections must be above a specified minimum S/N , where 
the S/N is given by equation (10). Including a threshold 
criterium, the number of MBHBs with S/N > s in an 



observation of duration r is then 



/"/ISCO /"/max Jt 

N T (> s) = R -iH s (Af) + R / 4-# s (A/ min ), 

Jfmin J J/„i. / 

(14) 

where 



In the second term of equation (14), which again accounts 
for the new binaries formed at the hardening radius, /max 
is the frequency reached after 3 years starting from / m i n , 
and the function H s (Af m - m ) is evaluated by integration 
of the S/N from / min to /. Given the exceedingly low 
value of / m in, this second term is totally negligible for an 
experiment such as LISA. 



3. HIERARCHICAL GROWTH OF MASSIVE BLACK HOLES 

The theory and method outlined in the previous sections 
allow us to fully characterize the expected contribution of 
MBHBs in the spiral-in phase to the LISA data stream, 
once the coalescence rate of MBHBs is specified. In this 
work a hierarchical structure formation scenario for the as- 
sembly and growth of MBHs in which seed holes form far 
up in the dark halo "merger tree" is assumed. We use ex- 
actly the same model discussed in Volonteri et al. (2003a, 
2003b) and in Paper I. Its main features are briefly sum- 
marized in this section. 

We track backwards the merger history of 220 parent ha- 
los with present-day masses in the range 10 11 — 10 15 M Q 
with a Monte Carlo algorithm based on the extended 
Press-Schechter formalism (see, e.g., Cole et al. 2000). 
Seed holes with m see( j = 150 M Q are placed within rare 
high-density regions (minihalos) above the cosmological 
Jeans and cooling mass at redshift 20. Their evolution 
and growth is followed through Monte Carlo realizations of 
the halo merger hierarchy combined with semi-analytical 
descriptions of the main dynamical processes, such as dy- 
namical friction against the dark matter background, the 
shrinking of MBHBs via three-body interactions, their co- 
alescence from the emission of gravitational waves, triple 
MBH interactions, and the effect of gravitational recoil. 
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Fig. 3. — Number of coalescences of MBHBs observed per year 
at z = per unit redshift. Our fiducial rate (thick solid line) is 
compared to a case in which the hardening timescale is increased 
by a factor of 3 (dotted line) or reduced by the same factor (dashed 
line) . 

Quasar activity is triggered during major mergers. We 
assume that the more massive hole accretes, at the Ed- 
dington rate, a gas mass fraction that scales with the fifth 
power of the host halo circular velocity (Ferrarese 2002). 

In a typical merger event, dynamical friction drives the 
satellite halo toward the centre of the new forming sys- 
tem, leading to the formation of a bound MBHB in the 
violently relaxed stellar core. As the binary separation de- 
cays, the effectiveness of dynamical friction slowly declines; 
the bound pair then hardens by capturing stars passing 
within a distance of the order of the binary semi-major 
axis and ejecting them at much higher velocities (gravi- 
tational slingshot). The heating of the surrounding stars 
by a decaying MBH pair creates a low-density core out 
of a preexisting stellar cusp, slowing down further binary 
hardening (see, e.g., Milosavljevic & Merritt 2001). If the 
hardening continues sufficiently far, GW emission takes 
over, driving the pair to coalescence. Figure 3 shows the 
number of MBHB coalescences per unit redshift per unit 
observed year predicted by our model: we expect ~ 60 co- 
alescences per year, the vast majority involving quite light 
binaries (Mi + M 2 < 10 5 M o ). The model was shown to 
reproduce rather well the observed luminosity function of 
optically-selected quasars in the redshift range 1 < z < 5 
and the evolution of the nuclear MBH mass density with 
cosmic time (Volonteri et al. 2003a), and to provide a 
quantitative explanation to the stellar cores observed to- 
day in bright ellipticals as a result of the cumulative erod- 
ing action of shrinking MBHBs (Volonteri et al. 2003b). 



Fig. 4. — Estimated LISA rms confusion noise (solid line), as 
the quadratic sum of the LISA instrumental single-arm Michelson 
noise hg (dotted line), the confusion noise from unresolved galactic 
(Nelemans et al. 200f , long-dashed line), and extragalactic (Farmer 
& Phinney 2003, short-dashed line) WD- WD binaries, and our es- 
timate of the confusion noise from unresolved MBHBs (thick-solid 
line) . 

4. NUMBER COUNTS 
4.1. Stochastic noise from MBHBs 

The customary definition of GW confusion noise level 
is the amplitude at which there is, on average, at least 
one source per frequency resolution bin. The frequency 
bin width is Af = 1/r, so the longer the observation, the 
smaller the noise. As pointed out by Cornish (2003), the 
crude "one bin rule" is much too simple to properly de- 
scribe a binary system. Using detailed information theory, 
Cornish (2003) shows that a GW background becomes un- 
resolvable when there is, on average, at least one source 
per eight bins. 

In the last decade a considerable effort has gone into 
quantifying the galactic and extragalactic confusion noise 
in the band 0.01 mHz - 1 Hz (e.g., Schneider et al. 2000; 
Freitag 2001; Nelemans et al. 2001; Farmer & Phinney 
2003). We have then applied the "eight bin rule" to asses 
the confusion noise associated with the evolving popula- 
tion of MBHBs, compared to the most recent estimates 
of the noises from galactic (Nelemeans et al. 2001, "one 
bin" rule, 1 year observation) and extragalactic (Farmer 
& Phinney 2003, "one-bin" rule, 3-year observation) white 
dwarf (WD) binaries. As shown in Figure 4, MBHBs pro- 
duce confusion noise at /<4 x 10~ 4 Hz. 

Figure 4 also shows the global LISA h Yms , along with 
separate contributions from different source populations. 
Though, as expected, MBHB stochastic noise dominates 
over WD- WD signals at low frequencies, it lies more than 
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an order of magnitude below the instrumental LISA sensi- 
tivity curve, and hence its contribution to the LISA h rms 
can be ignored. On the other hand, this hampers the possi- 
bility that LISA could take advantage of the MBHB noise 
to probe the cosmological evolution of such particular par- 
ent population. 



4.2. Mass function and redshift distribution 

We have divided the resolved sources into "merging" 
and "in-spiral" binaries (MBs and IBs, respectively). The 
former are those binaries that reach the ISCO during the 
duration of the observation with a signal above threshold. 




4 5 6 7 
log(M,) 



Fig. 5. — Mass distribution of the more massive member of MB- 
HBs resolved with S/N > 5 by LISA in a 3-year mission (solid 
line). The separate counts for MBs (short- dashed line) and IBs 
(long- dashed line) are also shown. 

These events are of particular importance, as they probe 
strong field effects and represent a unique chance of ob- 
serving the coalescence and ring-down phases of MBHBs. 
Resolved IBs, instead, do not allow a direct observation 
of the coalescence phase. These events arise from light 
binaries whose final coalescence phase lies below thresh- 
old, and from binaries of all masses with rigco > 3 
years at the very start of the observation. We expect 
very few IBs in this last stage anyway, because binaries 
above threshold have typically tisco < 3 years when 
the observation starts (see Fig. 2). An example of a 
resolved IB is represented by the lowest curve of Fig- 
ure 2. MBHBs in this class have an integrated S/N 
above threshold, though the coalescence phase occur at 
too high frequency to be directly observed by LISA. 




redshift 



Fig. 6. — Differential redshift distribution of MBHBs resolved 
with S/N > 5 by LISA in a 3-year mission. Line style as in Fig. 5. 

An obvious consequence of our classification is that MBs 
have larger mass and a lower redshift than IBs. The mass 
distribution of the most massive member of the binary 
Mi is shown in Figure 5. The differential and cumulative 
redshift distributions are plotted in Figure 6 and Figure 7, 
respectively. Detectable IBs are ~ 55 in total, and account 
for almost all of the MBHBs observed at z>7. Conversely, 
the rarer MBs (~ 33 in total) are confined to the red- 
shift interval 2<z<7. Figure 8 shows the average mass 
ratio Mi/M\ for MBs and IBs. As expected, this is larger 
for IBs, given their larger average redshift. The ratio de- 
creases at low redshift, as a consequence of the complicated 
merger history of host dark matter halos. The reason why 
the average mass ratio of MBs peaks at Mi/M\ ~ 0.15 
lies in the fact that both the probability of halo mergers 
(because of the steep P&S halo mass function) and the 
dynamical friction timescale increase with decreasing halo 
mass ratio. Hence, fast equal mass mergers are rare, while 
in more common unequal mass mergers it takes longer than 
an Hubble time to drag the satellite hole to the centre. 
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Fig. 7. — Cumulative (integral from z to oo) redshift distribution 
of MBHBs resolved with S/N > 5 by LISA in a 3-year mission. Line 
style as in Fig. 5. 



only MBs. The large GW-brightness of MBHBs is such 
that LISA will observe >90% of all coalescences occurring 
at z<5. The efficiency falls below 0.5 only for MBHBs 
at z>8. The efficiency to MBs only is, obviously, lower. 
Figure 9 shows that a space-based interferometer such as 
LISA can directly observe the final stage of the spiral-in 
phase of about half of all MBHBs coalescing at z ~ 5. 
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Fig. 8. — Mass ratio distribution of MBHBs resolved with S/N > 5 
by LISA in a 3-year mission. Line style as in Fig. 5. 

We can define the detection efficiency of a specific mis- 
sion as the number of observable events divided by the ex- 
pected number of coalescences in the same time interval. 
Figure 9 shows the global (MBs+IBs) detection efficiency 
for LISA and the efficiency considering as "detections" 



Fig. 9. — Detection efficiency (solid line), defined as the number 
of detected events (MBs+IBs) divided by the total number of coa- 
lescences in the same time interval, as a function of redshift. The 
efficiency considering only MBs as detections is also shown (dashed 
line) . 



5. DISCUSSION 

In this paper we have characterized the GW signal pro- 
duced by cosmological MBHBs, and we have then folded 
the signal into the LISA performance capabilities. We find 
that LISA should resolve more than 90% of all cosmolog- 
ical coalescences of MBHBs occurring at z<5. The detec- 
tion efficiency is already >0.5 for MBHBs at z ~ 8. We 
showed that the confusion noise from residual unresolved 
MBHBs is expected to be at least an order of magnitude 
below LISA instrumental noise. 

We have divided the resolved events into merging and 
in-spiral binaries. MBs are associated with systems at rel- 
atively low redshift involving heavy pairs (10 4 — 10 7 Mq). 
Their strong GW signal can be used to study the orbital 
evolution of the pair until the ISCO, allowing to test GR in 
extreme conditions. On the contrary, IBs are less massive 
pairs at higher redshift. Such systems can be generally 
observed, with moderate integrated S/N, only for a rela- 
tively short amount of time, from few weeks to few months, 
before the ISCO is reached. IBs are nevertheless impor- 
tant as they push the limit of observable MBHBs out to 
z ~ 12, and allow studies of the formation and assembly of 
seed holes of intermediate masses. Binaries at even earlier 
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epochs, while common in our model, are unfortunately too 
light to be observed by LISA with a relevant S/N. 

In our refined study, the <~ 20 "stationary sources" dis- 
cussed in Paper I, i.e. sources whose shift in frequencies 
A/ during the observing period is </, appear to be un- 
detectable. This is because a too optimistic sensitivity 
threshold for LISA was assumed in Paper I for frequen- 
cies below 0.1 mHz, where most of the stationary sources 
are expected. Indeed, an order of magnitude improvement 
of the LISA sensitivity below 0.1 mHz would lead to the 
detection of ~ 100 of such events. We stress here that 
our results may be sensitive to different model parameters 
and assumptions. As done in Paper I, we explore such 
possibility by running test models in which the hardening 
timcscalc th was divided and multiplyied by a factor of 3. 
The resulting coalescence rates are plotted in Figure 3. We 
find 50 coalescences per observed year in the "slow hard- 
ening case" and 78 in the "fast hardening" case, compared 
to 64 in our reference model. In terms of LISA number 
counts, the effects are small. In the slow hardening case, 
the coalescence rate decreases at z>9 as a large fraction of 
binaries has th longer than the then Hubble time tn- At 
lower redshifts it is th < tn anyway, the coalescence rate is 
basically unaffected with respect to the fiducial case, and 
so are the number counts (both for MBs and IBs). To 
obtain a significative reduction of observable sources, th 
must increase by a larger factor, so that also MBHBs in 
the range 5<z<10 will have th^tu- By increasing th by 
an order of magnitude, we find the number of coalescences 
per observed year decreasing to 30: in a 3-year observa- 
tion LISA would detect 20 (40) MBs (IBs). Note that, as 
increasing th results in a lower average redshift of coales- 
cences, this also increases the global detection efficiency. 
In other words, there are less sources, but a larger fraction 
of them is detectable. 

In the fast hardening case, more binaries can coalesce 
at early times, and the number of surviving MBHBs at 
6<z<12 decreases. This ultimately causes a slight reduc- 
tion in the number of IBs observable by LISA (49 in 3 
years, compared to 55). From Figure 3, it is also evident 
that, for z<5, the coalescence rate is almost identical to 
the standard case, implying that the number of detectable 
MBs will remain approximately the same. 

In Figure 8, we showed that IBs have an average mass 
ratio higher than MBs. To check whether this result de- 
pends on our assumption of equal-mass seed holes, we ran 
a test model with a flat initial mass function for the seed, 
10 < m S00( j < 500M Q . The resulting binary mass distri- 
bution relevant for LISA was basically unaffected. 

The vast majority of IBs are low mass systems at fairly 
high redshift. Their characteristic strain lies just above the 
LISA threshold at frequencies of the order of 10 -3 — 10 -2 
Hz (see Fig. 2) where confusion noises from unresolved 
galactic and extragalactic WD-WD binaries dominate the 



sensitivity curve (see Fig. 4). WD-WD confusion noise 
levels are difficult to compute because of the many un- 
certainties in stellar population synthesis models, and in 
estimating the fraction of binary stars in galaxies. In our 
fiducial sensitivity curve we have added to the LISA effec- 
tive noise the galactic WD-WD confusion noise computed 
by Nelemans et al. (2001), and the extragalactic WD-WD 
confusion noise estimated by Farmer & Phinney (2003). 
Note that Nelemans et al. (2001) assumed 1, rather than 
3, year integration, and both estimates employ the one bin 
rule. Using the eight bin rule we expect these noises to in- 
crease to some extent, but differences should be small. An 
alternative accurate estimation for galactic WD-WD con- 
fusion noise was performed by Hils & Bender (2000), who 
assumed a three bin rule. Using their noise level, which 
is somewhat higher than that computed by Nelemans et 
al. (2001), the number of IBs observed by LISA in 3 years 
slightly decreases, from 55 to 50. 

Compared to unresolved galactic WD-WD binaries, the 
uncertainties in the extragalactic WD-WD confusion noise 
are much larger, and may have a more relevant impact on 
number counts. The estimate of Schneider et al. (2000) 
lies nearly a factor of four above Farmer & Phinncy's "fidu- 
cial" model, and this leads to a significant reduction of 
observable IBs, from 55 to 44. The lowest estimate we 
could find in the literature was Farmer & Phinney's "pes- 
simistic" model, which is about a factor of four lower than 
their fiducial one; in this case we count 63 observable IBs. 

Another potential source of noise in the frequency range 
1-10 mHz is captures of compact objects (white dwarfs, 
neutron stars, and stellar-mass black holes) by MBHs in 
galaxy centers. Capture rates are quite uncertain, and 
estimates of relevant confusion noises span more than an 
order of magnitude in h c (see Barack & Cutler 2004 for 
a detailed discussion). We have estimated the impact of 
compact object captures on LISA number counts assum- 
ing the more optimistic rates calculated by Freitag (2001). 
The number of observable IBs decreases in this case to 43. 
More conservative rate estimates do not affect appreciably 
the counts. 

To summarise, in the context of our model, we can assign 
an approximate error of ~ 50% to the number of high-z 
MBHBs detectable by LISA. We conclude remarking that 
the bulk of detections involves binaries with masses in the 
interval 10 3 — 10 5 M Q , a range where black holes have never 
been observed. Genuine supermassive BH binaries, whose 
existence is more secure on observational grounds, appear 
too "heavy" for interferometers working in the band 0.01 
mHz - 1 Hz. 

We thank P. Bender and A. Vecchio for several enlight- 
ening discussions during the preparation of this paper. 
Support for this work was provided by NASA grant 
NNG04GK85G and by NSF grant AST-0205738 (P.M.). 
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